Analytic evaluation of the electronic self-energy in the GW approximation 

for two electrons on a sphere 
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The GW approximation for the electronic self-energy is an important tool for the quantitative 
prediction of excited states in solids, but its mathematical exploration is hampered by the fact that 
it must, in general, be evaluated numerically even for very simple systems. In this paper I describe 
a nontrivial model consisting of two electrons on the surface of a sphere, interacting with the normal 
long-range Coulomb potential, and show that the GW self-energy, in the absence of self-consistency, 
can in fact be derived completely analytically in this case. The resulting expression is subsequently 
used to analyze the convergence of the energy gap between the highest occupied and the lowest 
unoccupied quasiparticle orbital with respect to the total number of states included in the spectral 
summations. The asymptotic formula for the truncation error obtained in this way, whose dominant 
contribution is proportional to the cutoff energy to the power —3/2, may be adapted to extrapolate 
energy gaps in other systems. 

PACS numbers: 71.15.Qe, 71.45.Gm 



I. INTRODUCTION 

Accurate first-principles calculations of electronic ex- 
citations in solids are notoriously challenging, because 
the variational principle, which underlies ground-state 
schemes like Kohn-Sham density-functional theory^'^ or 
quantum Monte Carlo methods,^ cannot be exploited in 
this case. Therefore, many-body perturbation theory,* 
which is based on Green functions and allows an ex- 
plicit incorporation of relevant Coulomb-correlation ef- 
fects through a summation of the corresponding Feynman 
diagrams, is often the method of choice for quantitative 
ab initio descriptions of experimental spectroscopies. An 
especially fruitful realization of this framework for ac- 
tual electronic-structure calculations is the so-called GW 
approximation^ that yields quasiparticle band structures 
in much better agreement with photoemission data than 
density-functional theory with standard local or semilo- 
cal exchange-correlation functionals.^ From a mathemat- 
ical point of view, it constitutes an expansion of the exact 
nonlocal and frequency-dependent self-energy to first or- 
der in the dynamically screened Coulomb potential W, 
which describes the interaction between two quasiparti- 
cles formed by an electron or a hole together with its 
surrounding polarization cloud. The GW approximation 
is hence particularly suited for materials with weak to 
medium correlation strength, such as semiconductors or 
simple metals, but not for strongly correlated systems, 
where the quasiparticle picture breaks down. 

In spite of the undisputed success of the GW approxi- 
mation for the prediction of material properties and the 
interpretation of photoemission measurements, the de- 
bate over its best practical implementation shows no ten- 
dency of abating. One long-standing controversy cen- 
ters on the question whether the Green function used 
to construct the self-energy should be evaluated self- 
consistently or not.^"-'^^ Both approaches can be justified: 



The Green function obtained with full self-consistency is 
the variational solution that makes the Luttinger-Ward 
functional^^ for the total energy within the random-phase 
approximation stationary, while the more commonly 
applied non-self-consistent version follows naturally from 
the iterative solution of Hedin's coupled integral equa- 
tions for the self-energy with a mean-field treatment as 
the starting point. ^ 

Assessments of different variants of the GW approx- 
imation are typically based on a comparison of the nu- 
merical results for selected test systems. However, the 
adopted numerical procedures and limited computational 
resources inevitably necessitate additional simplifications 
whose impact on the final results is not always clear. For 
instance, the case studies investigating the effects of self- 
consistency for real solids reported so far have, in general, 
only used a restricted form of self-consistency that was 
limited to the quasiparticle energies in the denomina- 
tor of the Green function, the quasiparticle orbitals,^'^ 
or the diagonal part of the Green function. Even at 
the level of the standard non-self-consistent GW approxi- 
mation, far-reaching additional simplifications are widely 
employed. Among these are the pseudopotential approx- 
imation, which leads to small but systematic deviations 
from all-electron results due to the inexact core-valence 
partitioning and the use of pseudo wave functions, ^""^^ as 
well as plasmon-pole models or other simplified screening 
functions.^ A range of less obvious but equally important 
factors like the proper treatment of anisotropic screening 
in noncubic systems^"'^^'^^ or problems resulting from an 
incomplete basis set in all-electron calculations^^'^'* have 
also been emphasized. 

One point that has increasingly come into focus in this 
context is the convergence behavior with respect to the 
number of unoccupied states included in the construction 
of the Green function G and the screened Coulomb inter- 
action W. Much attention was raised after an early all- 
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electron calculation,"'^^ which claimed much smaller semi- 
conductor band gaps than established pseudopotential 
results, was put into question due to alleged incomplete 
convergence with respect to this parameter, leading 
to further detailed studies^^'^^'^^ that illustrated the slow 
convergence rate of the quasiparticle energies for a wide 
variety of materials. As an extreme example, Shih et al?"^ 
reported that bulk zinc oxide required thousands of unoc- 
cupied bands to achieve satisfactory accuracy in a calcu- 
lation based on pseudopotentials and a plane- wave basis 
set. It was subsequently argued that this peculiar behav- 
ior was caused by the particular choice of plasmon-pole 
model used in Ref. 32 and that far fewer bands are in fact 
required if the full frequency-dependent screening func- 
tion is properly constructed within the random-phase 
approximation.^^ However, a similarly slow convergence 
was again observed in an all-electron calculation for zinc 
oxide that not only avoided plasmon-pole models but also 
the additional pseudopotential approximation.^^ Paral- 
lel to these developments, different approaches were pro- 
posed to circumvent or at least alleviate the convergence 
problem. These include the replacement of all high-lying 
empty states by plane waves, the extrapolar method, "^^ 
in which merely a small number of unocciipicd states are 
treated explicitly and a common energy denominator is 
assigned to the remainder, so that the closure relation 
can be applied, the Lanczos-chain algorithm, as well as 
the effective-energy technique^^ and methods based on 
the self-consistent Sternheimer equation, which are for- 
mally exact despite only involving occupied states. While 
these novel schemes undoubtedly hold great potential, 
practical applications are not yet widespread, in part be- 
cause not all popular computer codes support them at 
present. As a consequence, the majority of GW calcula- 
tions still rely on traditional procedures and suffer from 
the problem of slow convergence. 

All of the above issues are relevant on the energy scale 
of several tenths of an electron volt that matters for 
the comparison between different implementations of the 
GW approximation and with experiments, but their con- 
trol is difficult in practice due to complex interdependen- 
cies. Therefore, model systems that permit numerically 
exact or, ideally, analytic solutions play an important 
role for developing and testing approximation schemes 
within many-body perturbation theory, but even the ho- 
mogeneous electron g frequently employed model in 
solid-state physics, can only be treated numerically in the 
GW approximation. Furthermore, with no experimental 
measurements or independent theoretical benchmark re- 
sults, even the basic question whether the true occupied 
band width in the range of metallic densities is smaller 
than that of free electrons, as predicted by the standard 
non-self-consistent GW approximation, ^'^^ or larger, as 
obtained when full self-consistency is included,^'* is not 
yet finally settled. Calculations that go beyond the GW 
approximation and attempt to incorporate the combined 
effects of self-consistency and vertex corrections remain 
inconclusive, because the results depend on the choice of 



vertex function and details of the implementation.'''^^''*'' 
In this situation, more tractable few-electron systems are 
of considerable interest. 

The first nontrivial system for which the self-energy at 
the GW level can be derived analytically was a Hubbard 
model with four sites in a tetrahedral arrangement and 
two electrons, originally used as a counterexample to 
demonstrate the violation of particlc-mimbcr conserva- 
tion in the non-self-consistent GW approximation before 
a more general investigation of this problem based on 
symmetry arguments.'*^ The analytic solvability was im- 
portant in this case, because it proved unequivocally that 
the quantitative deviation was genuine and not due to nu- 
merical inaccuracies. A related but even simpler two-site 
model with a pair of electrons can be treated analytically 
in the same way.**^ Lattice models with a wider range of 
parameters, for which the GW self-energy is accurately 
obtainable by numerical means, were also employed in 
several studies. *°'^'*"'*^ The properties of Hubbard mod- 
els deviate in many respects from those of real materials, 
however, and conclusions from such comparisons cannot 
always be directly transferred to the ah initio realm. 
Most importantly, the local on-site interaction differs sig- 
nificantly from the actual Coulomb potential and leads 
to a dominance of short-range correlation effects. These 
are not well described by the GW approximation, which 
mainly accounts for the long-range screening of charge 
carriers. Furthermore, the restricted Hilbert space does 
not allow us to address problems like the convergence 
behavior with respect to the number of empty states. "^^ 
Peculiar symmetries, such as that between occupied and 
unoccupied states in the two-site model at half filling,*^ 
which are not obeyed by real solids, may also have an 
influence on the results. For completeness, it should be 
mentioned that the polaron model of individual electrons 
coupled to an external boson field can also be treated 
analytically,^*^ but its usefulness as a test system for the 
GW approximation is even more limited, as there is no 
explicit renormalizable electron-electron interaction. 

For future methodological investigations I here propose 
a better suited continuum system consisting of two elec- 
trons confined to the surface of a sphere, and I show that 
the self-energy within the standard non-self-consistent 
GW approximation can be derived entirely analytically. 
In contrast to the previously considered Hubbard mod- 
els with the same property, '''^'''^ the electrons interact 
with the normal long-range Coulomb potential, and there 
is an infinite Hilbert space of single-particle wave func- 
tions whose eigenvalues are not bounded from above. 
Therefore, the performance of particular approximation 
schemes should be more indicative of applications to real 
systems. The model considered here can be regarded as 
a two-dimensional homogeneous electron gas in a closed 
curved space, whose density depends on the radius of 
the sphere. As in the three-dimensional electron gas, 
the correlation is weak at high densities (small radius) 
and becomes strong at low densities (large radius),''^ so 
that different regimes can be explored within the same 
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framework. The system is also quasi-exactly solvable, 
which means that some exact eigenvalues, although not 
the complete spectrum, are known analytically."'''' 

This paper is organized as follows: In Sec. II the system 
is mathematically defined and discussed in more detail, 
before the analytic expression for the self-energy in the 
non-self-consistent GW approximation is derived in Sec. 
III. Then in Sec. IV the convergence behavior with re- 
spect to the number of empty states is analyzed, leading 
to an analytic formula describing the asymptotic depen- 
dence on the cutoff energy. The conclusions are summa- 
rized in Sec. V. Unless otherwise noted, Hartree atomic 
units are used throughout. 



II. MODEL DESCRIPTION 

The system considered here consists of two electrons 
on the two-dimensional surface of a sphere with radius 
R. Their positions are expressed in spherical coordinates 
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in terms of the polar angle and the azimuthal angle <p. 
A homogeneous positive surface charge density 2/{A-kB?), 
which gives rise to the attractive electrostatic potential 
— 2/i?, is included to ensure overall charge neutrality. If 
the self-interaction of the positive charge background is 
also taken into account, then the Hamiltonian becomes 
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which is related to the angular-momentum operator L. 
Its eigenfunctions are the spherical harmonics 
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The eigenstates of the Hamiltonian (2) must in general 
be determined numerically, but for certain discrete radii 
individual analytic solutions are known. ^° For example, 
for the particular value R = a/3/2, which corresponds to 
intermediate correlation strength, the exact ground-state 
energy for two electrons is Eo(2) = 1 — 4/-\/3. 

Due to symmetry requirements, the ground-state elec- 
tron density is evenly distributed on the spherical surface. 
In density-functional theory the one-particle Hamiltonian 
of the auxiliary Kohn-Sham system thus takes the form 
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with a constant exchange-correlation potential V^c- The 
Hartree potential is exactly canceled by the electrostatic 



potential of the positive charge background. Prom Eq. 
(4) the Kohn-Sham orbitals are 
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where the normalization is chosen with respect to the 
two-dimensional integral over the spherical surface 
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and the corresponding eigenvalues 
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are independent of the magnetic quantum number to. In 
the ground state the lowest spin-degenerate Kohn-Sham 
orbital is doubly occupied, whereas all others are unoc- 
cupied. As the eigenvalue of the highest occupied orbital 
also equals the negative of the ionization potential in ex- 
act density-functional theory,^^ the exchange-correlation 
potential can be determined from the difference 
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between the ground-state total energy £"0(2) of the true 
interacting two-electron system and the energy Eq(1) — 
of the corresponding ionized one-electron system, if the 
former is known. Thus for R = y/3/2 the exact exchange- 
correlation potential is V^c = 1 — 4/\/3. Alternatively, for 
free electrons V^^ is set to zero. Analogous to the three- 
dimensional homogeneous electron gas, these two choices 
differ only by a trivial energy shift. 



III. DERIVATION OF THE SELF-ENERGY 

The first ingredient required for the construction of the 
GW self-energy is the Kohn-Sham Green function. In the 
following, all functions will be projected onto the orbitals 
(6), which form a complete set. The Green function is 
diagonal in this basis, and the diagonal matrix elements 
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are furthermore independent of to. The symbol rj always 
denotes a positive infinitesimal. Since the wave function 
2/00 = (47rii2)~^/2 of the only occupied Kohn-Sham state 
is constant, the polarization function becomes 
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with the definition = — eo and a factor 2 for the spin 
summation. The representation of the Coulomb potential 
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follows from the addition theorem for the spherical har- 
monics. It is immediately clear from this formula that the 
interaction strength grows with R while the level spacing 
between the eigenvalues (8) is simultaneously reduced, so 
that the system becomes more and more strongly corre- 
lated with increasing sphere radius. As the polarization 
function and the Coulomb potential are both diagonal in 
the chosen basis, the matrix elements of the dynamically 
screened interaction 



with the correlation part 
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can be calculated by means of a simple scalar renormal- 
ization. The poles are located at 
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The exchange part of the self-energy, which is frequency 
independent and equals the nonlocal exchange potential 
in Hartree-Fock theory, is obtained as 
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while the correlation part of the self-energy is given by 
a convolution of the Green function and the screened 
interaction. The matrix elements are given by 
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These overlap integrals are in fact real valued and zero 
unless mi + m2 = m. They are written conveniently in 
terms of the Wigner 3-j symbols^'^ as 
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FIG. 1. Right-hand side -I- -I- E^(a;) - 14c of the quasi- 
particle equation (25) for the occupied state (£ = 0, above) 
and the first unoccupied state (I — 1, below) for R — \/3/2. 
The solutions of the quasiparticle equation correspond to the 
intersections with the diagonal (dashed line). In addition 
to the actual quasiparticle, marked by a circle, there are an 
infinite number of satellite resonances for each value of I. 



It then follows from the orthogonality relation 
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that the correlation part of the self-energy, like all other 
quantities considered in this section, is diagonal in the 
Kohn-Sham basis and that the diagonal elements 

(2^i-M)(2£2 + l) f £1 £2 £^^ 
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are again independent of m. Carrying out the remaining 
contour integration eventually yields 
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This general formula can now be exploited to derive the 
self-energy correction of individual orbitals. Of particular 
interest are the highest (and only) occupied state 
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as well as the lowest unoccupied state 
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Finally, the energy spectrum of electronic excitations is 
given by the solutions of the quasiparticle equation 
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In Fig. 1 both sides of this nonlinear equation are shown 
for ^ = (above) and i = 1 (below) for a sphere radius 
of i? = VS/2. The possible excitations correspond to the 
intersections of the two curves. Evidently, there are an 
infinite number of satellite resonances for each quantum 
number £, stemming from the poles of the self-energy, 
in addition to the principal quasiparticle state, which is 
marked by a circle. It should be noted that the poles 
are well separated at the positions ee + z^ with integer £ 
for the occupied state, while the self-energy of the unoc- 
cupied state features pairs of poles located very close to 
each other at + ze-i and e^-i + Zi. These stem from 
the second and third term on the right-hand side of Eq. 
(24), but the splitting cannot be properly resolved on the 
scale of the figure. As a rule, the multiplicity of the satel- 
lites splitting increases with the quantum number of the 
quasiparticle state. If the self-energy is linearized around 
the Kohn-Sham eigenvalues e^, as is common in ab initio 
calculations, then the quasiparticle energies are given by 
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where the renormalization factors 
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(27) 



specify the weight of the quasiparticle resonance in the 
spectral function. For R = •\/3/2 the values are Zq « 0.94 
and Zi w 0.95. 

The energy gap Egap = e^''— ej}^ thus obtained is shown 
in Fig. 2 as a function of the sphere radius R relative to 
the Kohn-Sham eigenvalue gap -E-^p = ei — eq. In addi- 
tion, the exact value, which is defined as the difference 
I — A between the ionization potential I = Eq{1) — Eo{2) 
and the electron affinity A = Eo{2) — Eq{3), is displayed 
for comparison. The ground-state energies Eo(2) for two 
electrons and Eq (3) for three electrons are obtained from 
a numerical diagonalization with a full set of Slater deter- 
minants constructed from the orbitals (6), converged at 
the scale of the figure. For small radii, i.e., high electron 
densities, the dynamics of the system are dominated by 
the kinetic energy, and the quasiparticle energy gap ap- 
proaches that of nonintcracting electrons. With increas- 
ing sphere radius, correlation effects become stronger, 
and the quasiparticle energy gap widens relative to the 
eigenvalue gap of the nonintcracting Kohn-Sham system. 
Up to an intermediate correlation strength of i? ~ 1 , the 
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FIG. 2. Energy gap -Egap = e^'' — (-q'^ between the occupied 
and the lowest unoccupied quasiparticle state as a function 
of the sphere radius R relative to the Kohn-Sham eigenvalue 



gap Eg 



£() (dashed line). The exact numerical value, 



obtained from the difference between the ionization potential 
and the electron affinity, is also shown (solid line). Small R 
correspond to weak and large R to strong correlation. 



GW approximation is in excellent quantitative agreement 

with the exact numerical value, but it underestimates the 
further rapid increase of the gap in the strong-correlation 
regime at larger radii. Altogether, the behavior of the 
GW approximation for this model system hence accords 
completely with its performance for real solids. 



IV. ASYMPTOTIC CONVERGENCE 

The analytic expression (22) for the self-energy derived 
above includes a double infinite summation over the an- 
gular quantum numbers £i and £2, reflecting the spectral 
sums over unoccupied eigenstates in the Green function 
and the screened interaction, respectively. In practice, 
such sums must be truncated at a certain cutoff energy, 
so that quantitative deviations from the true results are 
incurred as a consequence. In the following I analyze the 
asymptotic convergence with respect to the number of 
unoccupied states for this system. 

Owing to the truncation, the exact matrix element of 
the correlation part of the self-energy in Eq. (26) is re- 
placed by Y^Kei), which includes only those terms of Eq. 
(22) where both £1 and £2 arc smaller than a particular 
finite quantum number Lcut • If the difference is denoted 
by Ai = S^(e£) — ^^(ef), then the associated error in the 
quasiparticle energies due to the truncation equals Z^Af. 
In principle, there is also another distinct error that stems 
from the approximate evaluation of the renormalization 
factors (27), which are in practice obtained from the trun- 
cated S^(u;) instead of the exact S^(w), but it turns out 
that the resulting additional deviation is proportional to 
the cutoff energy to the power —2 and hence not relevant 
for the following analysis of the leading-order corrections. 
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In accordance with Eq. (23), the truncation changes 
the self-energy matrix element for the occupied state by 
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If the addends on the right-hand side are expanded in 
inverse powers of € according to 
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then the entire expression can be rewritten as a sum of 
Hurwitz zeta functions 
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with positive integer exponents n. From the asymptotic 
behavior of the Hurwitz zeta function^^ 
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In the last step a cutoff energy 
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which is half-way between the energies of the states with 
-^cut — 1 and Lent , was inserted. An analogous calculation 
for the first excited state based on Eq. (24) yields 
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As the terms proportional to L~^^ fail to fortuitously can- 
cel in this case, the final expression retains a nonvanishing 
contribution with the cutoff energy to the power —3/2. 
The truncation error of the energy gap between the occu- 
pied and the lowest unoccupied quasiparticle state hence 
exhibits the asymptotic behavior 
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The leading term is proportional to E~^^ but of small 
absolute magnitude, because the renormalization factors 
(27) in the vicinity of the fundamental gap show very 
little variation. In fact, if the quasiparticle energies are 
evaluated with Z^ ~ 1, as has been advocated by some 
authors, 2® then this term vanishes exactly. In practice, 
the convergence of the gap is hence dominated by the 
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term proportional to -Eput ■ This remains true for other 
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FIG. 3. Convergence of the matrix elements of the correlation 
part of the self-energy E^(ef) for the occupied state [l — 0) 
and the first unoccupied state {£ — 1) as well as their contri- 
bution ZiS5(£i) ~ ■^oSo(eo) to the energy gap for R — y/3/2. 
The dashed lines indicate the asymptotic behavior given by 
the analytic formulas (32), (34), and (35). In the latter case, 
only the term proportional to E^^^ is included. 



transitions; for example, between the occupied state and 
the second unoccupied state, whose truncation error 
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can be derived along the same lines. If the cutoff energy 
is not chosen as in Eq. (33) in the center but elsewhere 
in the interval between the states with quantum num- 
bers Lcut — 1 and Lent, then the formulas for all Ai are 
modified by an additional identical term proportional to 



E^ 



-3/2 



As a visualization. Fig. 3 displays the matrix elements 
^jiee) for the occupied state {£ = 0) and for the first 
unoccupied state {£ = 1) as a function of the cutoff en- 
ergy for R = VS/2. The steplike variation refiects the 
discrete nature of the Kohn-Sham eigenvalue spectrum. 
The dashed lines indicate the asymptotic expansions 
1^1{ee) + Ai{Ecnt), where the analytic expressions for Aq 
and Ai from Eqs. (32) and (34) with terms up to the or- 
der i?cut^^ ^I'S used. The difference Zi'S'Kei) — ZoS§(eo), 
which equals the contribution of the correlation part of 
the self-energy to the gap, is also shown together with the 
asymptotic formula Zii;^[(ei)-ZoSo(£o) + Agap(i?cut)- In 

3/2 

this case, only the term proportional to iS^ut '^^ ^'i- (35) 
is considered, because the difference Zi — Zq ~ 0.01 is 
so small that the leading term is almost completely sup- 
pressed; the change resulting from its inclusion would not 
be discernible on the scale of the figure. The asymptotic 
convergence of the energy gap is evidently well described 
in this way. 

In actual ab initio calculations for real solids, the high 
computational cost often precludes a convergence to the 
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desired accuracy. The contribution of the omitted high- 
lying unoccupied states may be approximately included 
within the extrapolar method^"'^ or the effective-energy 
technique,^'' but straightforward extrapolation to arrive 
at the limiting values would seem the most natural course 
of action where implementations of these schemes are not 
available. In practice, however, a direct extrapolation of 
the quasiparticle energies has only been attempted very 
rarely due to uncertainties about the proper asymptotic 
formula. In Ref. 34 an expression of the form 

where N is the number of bands and a, b, and iVo are fit- 
ting paramciters. was ciniploycd ad hoc to extrapolate the 
limiting value b for the band gap of zinc oxide, whereas 
fitting functions with different powers of N have been 
used in other studies. As the number of bands increases 

3/2 

proportional to E^l^^ at high cutoff energies, this work 

suggests that the form f{N) — & ~ A''"-^ ~ E^ut^ guessed 
in Ref. 34 is indeed correct. The suppression of the lead- 
ing order in Eq. (35) also provides an explanation why 
band gaps are often observed to converge faster than the 
individual quasiparticle energies. 

V. CONCLUSIONS 

In this paper I have described a nontrivial model sys- 
tem of two interacting electrons on a sphere for which 
the self-energy in the GW approximation without self- 
consistency can be evaluated analytically. This corre- 
sponds to the standard approach taken in virtually all 
actual ab initio calculations. As the relevant character- 
istics, such as the long-range Coulomb potential and the 
infinite Hilbert space, are the same as in real materials, 
the system appears better suited to explore the proper- 
ties of the GW approximation than previously employed 
analytically solvable lattice models. Indeed, the results 
presented here demonstrate an analogous performance 
as for real materials: The GW approximation corrects 
the underestimation of the fundamental energy gap in 
Kohn-Sham density-functional theory and yields accu- 
rate quantitative results for low to intcirmediate correla- 
tion strength, but fails in the strong-correlation regime. 
For this reason, it suggests itself natural testing 

ground to study extensions beyond the standard GW ap- 
proximation that are designed to describe the self-energy 



of strongly correlated systems. Even if a purely analytic 
treatment is then no longer possible, the computational 
cost will be much smaller than for real materials, allowing 
a highly accurate evaluation without the apparent artifi- 
cialities and the parameter dependence of typical lattice 
models, which have repeatedly been chosen to study the 
influence of vertex corrections and self-consistency in the 

pg^g^ 10,11,43,46,56,57 

Here this systcim was used to study the convergence 
of the self-energy with respect to the number of empty 
states included in the spectral summations. The results 
not only confirm previous empirical observations that 
transition energies converge faster than individual quasi- 
particle states due to a partial error cancellation, but the 
asymptotic expansion also demonstrates that the gap be- 
tween the highest occupied and the lowest unoccupied 
state approaches its limiting value with an error that is, 
for practical purposes, proportional to the cutoff energy 
to the power —3/2. Although a more general study of 
the asymptotic behavior is highly desirable, there is no 
indication that the dependence on the cutoff energy ob- 
tained here is due to specific details of this model. In- 
deed, the truncation error of the quasiparticle band gap 
in solids appears to exhibit the same exponent. "^^ If con- 
firmed, this would enable practical direct extrapolation 
schemes with an appropriate fitting function, implying 
enormous potential benefits for computationally expen- 
sive ab initio calculations. Furthermore, if some way was 
known to determine the relevant prefactor based on gen- 
eral characteristics of the material in question, then even 
an a posteriori correction without the ambiguities of nu- 
merical fitting procedures would be possible. Such an a 
posteriori scheme exists, for example, to extrapolate the 
self-energy correction of the band gap from the repeated- 
slab approximation with finite vacuum buffers and three- 
dimensional periodicity to the limit of an isolated slab, 
greatly accelerating the convergence of GW calculations 
for thin films with respect to the supercell size.^® In this 
sense, it is hoped that the present paper sparks further 
fruitful work along the same lines. 
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